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I.  INTRODUCTION 


Various  approaches  have  been  taken  to  analyze  theoretically  the 

performance  of  cw  diffusion -type  HF/DF  chemical  lasers.  These  include  the 

development  of  computer  programs  that  use  boundary  layer  approximations 

1-  3 

in  the  mixing  region,  an  "exact"  analytical  flame- sheet  solution  for 

4 . .5 

laminar  mixing  in  the  limit  of  zero  power  or  optical  saturation,  and 

solutions  in  which  simplified  mixing  and  chemical  kinetics  models  are 
6 7 

used.  ’ The  latter  do  not  satisfy  the  equations  of  motion,  even  in  a limiting 
sense,  but  do  describe  the  basic  behavior  of  the  cw  diffusion- type  chemical 
laser  in  terms  of  relevant  normalized  parameters. 

A flame-sheet  diffusion  model  was  used  by  Mirels,  Hofland,  and  King.  ^ 
Laminar  and  turbulent  mixing  were  analyzed  as  well  as  oscillator  and 
amplifier  configurations.  In  order  to  obtain  simple  closed-form  solutions, 

a two-level  vibrational  model  was  used. 

7 8 

Broadwell  extended  the  analysis  of  Emanuel  and  Whittier  to  include 

effects  of  diffusion  on  HF  laser  performance.  A scheduled-mixing  model  was 

used;  turbulent  mixing  and  an  oscillator  optical  configuration  were  assumed; 

multiple  vibrational  levels  were  considered. 

In  the  present  study,  the  analysis  of  Mirels,  Hofland,  and  King^  is 

extended  to  include  multiple  vibrational  levels.  Arbitrary  mixing  rates  and 

chemical  rate  coefficients  are  considered.  Solutions  for  zero  power  and 

amplifier  cases, as  well  as  expressions  for  total  oscillator  output  power, 

6 7 

are  obtained.  In  addition,  the  flame- sheet  and  the  scheduled-mixing  models 


r 


are  compared.  The  formulation  presented  herein  can  be  used  for  any 

diffusion- type  chemical  laser.  Both  cold- reaction,  e.g.  , + F -•  HF(v)  + H), 

and  chain-reaction,  e.g.  , + F - HF(v)  + H,  H + F^  - HF(v)  + F)  HF 

lasers  are  considered.  These  developments  go  beyond  the  treatment  of 
7 

Broadwell  in  that  he  assumed  a cold  reaction  HF  laser,  turbulent  mixing, 
and  particular  numerical  values  for  chemical  rate  coefficients. 
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II.  THEORY 


In  the  following  sections,  the  rate  of  formation  and  deactivation  of  excited 
species  in  a cw  diffus ion-type  chemical  laser  is  discussed.  Equations  that 
describe  the  performance  of  amplifier  or  oscillator  configurations  are  then 
noted.  In  order  to  fix  the  ideas,  a cold-reaction  HF(DF)  laser  is  considered  in 
Sections  II  and  III.  The  extension  to  a chain-reaction  HF(DF)  laser  is  noted  in 
Appendix  A. 

A.  BASIC  EQUATIONS 

The  mixing  model  of  Mirels,  Hofland,  and  King^  is  illustrated  in  Fig.  1. 
The  flow  downstream  of  a single  oxidizer  (F)  semichannel  is  considered.  The 
reactants  are  assumed  to  be  premixed,  but  do  not  start  to  react  until  a specified 
flame  sheet  location  is  reached.  The  effect  of  diffusion  is  accounted  for  by  the 
flame  sheet  shape.  The  flame  sheet  ordinates  (x^,  y^)  are  given  by 

yf  = yf("o>  ^o-^D 


= w 


> X 


D 


(lb) 


where  w is  the  width  of  the  semichannel,  and  x^^  is  the  axial  distance  required 
for  the  flame  sheet  to  reach  the  channel  centerline.  Flow  velocity,  pressure, 
and  temperature  are  assumed  to  remain  constant.  The  semichannel  is  assumed 
to  have  a unit  height. 


Figure  1.  Mixing  Model  of  Mirels,  Hofleind,  and  King. 

Reactants  are  premixed.  Reaction  starts  at  Xq 
i.  e. , at  intercept  £f  strearn_tube  with  flame 
sheet  y£  = yf(xo).  = I + + I^.  (a)  Typical 

physical  flow;  (b)  Corresponding  model  semi- 
channel flow. 


Let  denote  the  number  density  (mcde/cm  ) of  molecules  at  vibrational 
level  V.  Vibrational  levels  from  v = 0 to  v = v^  are  considered.  The  molecules 
in  each  vibrational  level  are  assumed  to  be  in  rotational  equilibrium  at  a 
temperature  equal  to  the  local  translational  temperature  T.  Only  P -branch 
transitions  (v  + 1,  J - 1 -»  v,  J)  are  considered.  The  gain  associated  with 
these  transitions  can  be  expressed  as 


g T - i « - A T-n  ) 

v,  J V,  J v+ 1 V,  J V 


where 


A T = e 
V,  .T 


■2JT^/T 


The  value  of  a at  line  center  of  a Doppler -broadened  P-branch  transition  can 

V,  J 


be  found  from^ 


47  1/2  3/2  2 -J(J-1)T_/T 

CT  - = 3.29  X "^T  ^Je  ^ 

J K V,  J 


wlic’i-e,  for  an  HF  lasing  moh-cule,  a 1,  I’j^  30.  16  K,  M 20  g/mole, 


2 3 8 3 

M - X 10-^  = 0. 96  [1  + 0.063J]  erg-cm-^ 


^/M^  ^ = 1 + V - 0. 01  V- 
V,  J o,  J 


*^Q  1 ~ 5.66  X 10  (400/T)^ ^^cm^/mole 
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Equations  (2d)  and  (2o)  are  correlations  of  estimates  provided  by  Herbelin  that 
are  believed  to  be  accurate  to  within  about  10%  for  I < J <16  and  v < 6. 


For  convenience,  the  symbols  r a ,,  and  A . are  now  denoted  g , 

V,  J V,  J V,  J 

Oy*  and  Ay,  respectively;  i.  e.  , the  dependence  on  J is  not  displayed.  It  is 
assumed  that,  when  lasing  occurs,  only  one  transition,  with  intensity  denoted 
J = is  permitted  between  v + 1 and  v.  In  order  to  obtain  closed- 
form  solutions,  vibration-vibration  (v-v)  collisions  are  n<‘glected.  The  /ate 
of  change  of  Oy  in  each  stream  tube  can  then  be  expressed  as 

^[(^)v  -(t)v-i1  C3) 


where  u(dny/dx)p  is  the  rate  (mole/cm  /sec)  of  formation  of  Oy  by  the  chemical 

"pumping"  reaction,  as  discussed  in  Appendix  A.  The  second  term  on  the  right 

side  of  Eq . (3)  is  the  net  rate  of  formation  of  Oy  by  collisional  deactivation 

(see  Appendix  A)  and  by  spontaneous  emission.  The  quantities  and  are 
^ cd  sp 

the  rates  (sec  ) of  collisional  deactivation  and  spontaneous  emission, 

respectively,  from  v to  v- 1 levels.  Note  that  k^°^  = k^“^  = 0.  Generally 

sp  cd 

(v) / ( v) 

^'sp'^^cd  ^ lasers.  The  last  term  on  the  right  side  of  Eq.  (3)  is  the 

net  rate  of  formation  of  Oy  due  to  stimulated  emission. 
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We  now  introduce  the  following  normalized  variables 


C = xk^‘ >/u 


- 1 


1 = a I 

V V VI 


N.  = (n  w)  / n.dy 
1 r ' ’ 


V I cci  sp  1/  cd 


G = (a 

V V 


JkA'l) 

,->r' 

"r'"'  / 


N , , - A N 

v+ 1 V V 


(4) 


The  quantity  n^  is  a reference  number  density,  which,  for  a cold  - reaction 
HF  laser,  is  taken  to  equal  (np)„-  Equation  (3)  is  integrated  with  respect  to  y, 
between  the  limits  y = 0 and  y = y^.  If  it  is  assumed  that  r^  is  a constant,  I^ 
is  a function  only  of  x,  i.e.  , a mean  value  is  used  at  each  station,  and  n^  = 0 
at  the  flame  sheet,  the  result  is 

dN  /dC  + (r  + I , + A I )N  - (r  ^,  + 1 )N  - (AIN)  , 

V '*  V V-1  V V V V+1  V v+1  v-1 

= (dN^/dOp  (5a) 

witli  tlie  boundary  condition  N (0)  - 0.  An  alternative  form  of  Eq.  (5a) 
that  retains  G is 

V 

dN  /dr  - (r  ^,N  - r N ) - (I  G - I ,G  = (dN  ! dr\  (5b) 

V ^ v+1  v+1  V V V V v-1  v-1)  V 
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Let  Since  can  be  changed  only  by  the 

v=0 

it  follows  that 


(dN^/dOp 


For  a cold -reaction  HF  laser,  Eq.  (A-5), 

(dN  /dr)  = a K,N^ 
V ^'p  V 1 F 


where 


K^N^  = K 


^ r> 


K. 


C dY 

^o o 

dC„ 


dC, 


= dY/dC  (Kj-^«>) 

Here  Y^  = y^(x^)/w,Y  = y^(x)/w,  \ = Z ^ 

^ v=0  ^ P P 

n^  = (np)^.  The  rate  of  change  of  the  total  population  is, 
and  (6  a) 


= K,Nj, 


pumping  process, 

(5c) 

{6a) 

(6b) 

(6c) 

= k /k^  and 
1 p cd 

from  Eqs.  (5c) 

(7) 


Conservation  of  F atoms  yields 


’I 


It  follows  that 


’^ToT 


^iCo 


7 


e Y dC 
o ® o 


(K^— -co) 


(9a) 

(9b) 


The  flame  sheet  shape  can  be  expressed  in  the  form 


Y = (C/Cq) 


Y = 1 


m 


Y = 1 
o 


m 


(C  s Cj^) 


(C  > Cq) 


(10a) 


(10b) 


where  m = 1 /2  and  1 for  laminar  and  turbulent  mixing,  respectively. 

Now  let  denote  the  number  of  photons  emitted  up  to  station  f , per 
initial  n^  particle,  from  the  transition  v + 1,  J - 1 — ► v,  J.  Thus 


= (un 


/-X  y^(x) 

^w)'  / dx  / ( 

•'n  ^ n 


(gl /f )^dy 


(11) 


i: 


1 G dC 

V V 


:i2) 
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▲ 


The  quantity  is  a n^easure  of  the  local  emitted  photon  flux  per 

initial  n^  particle.  Let  P denote  the  net  output  power,  per  semichannel,  up 
to  station  Q . It  follows  that 


P 


v^-1 

un  w ^ 0 

r V 

v=0 


(13) 


Vf-1 

When  a mean  photon  energy  € is  used,  P = un  w60,  where  0 = y'  0 is  the 

r V 

v=0 

net  photon  output  per  initial  n^  particle  up  to  station 
B.  ZERO-POWER  AND  AMPLIFIER  CASES 

When  is  specified,  Eq.  (5a)  is  linear  and  generally  can  be  integrated 
in  closed  form  in  order  to  find  as  a function  of  Q . The  net  gain  per  semi- 
channel is  then  found  from  G = N , - (AN)  . 

For  the  zero-power  case  (I^  = 0),  Eq.  (5a)  yields 


r C 

N e = a K, 
V V 1 


( 


^ c 

\rV 


f 

J o 


r C 


(14) 


Equation  (14)  is  integrated  successively,  starting  with  v = v^.  Note  that 
r^  = 0,  Expressions  for  N^  and  are  given  in  Appendix  B. 
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r 

j 

It  is  of  interest  to  expand  Eq.  (5a)  in  a power  series  about  Q = 0. 
Consider  the  case  Kj  — ►«,  Y = (^/C  found  that 

N/Y  a + C(m+t)''|(>-,,+I)a,,+  i 

V V vlt  V v)l  V-1 

- (r  + 1 , + A I )a  ] - _ + O(C^)  (15a) 

V V- 1 V V V C =0  ^ 

For  the  zero-power  case 

N^/Y  = a^  + C(m  + l)‘‘l(ra)^^j  - (ra)J  + 0(C  ^)  (15b) 

The  corresponding  gain  per  seniichannel  is 

= K+l  - ^vV  + 

- A^|(ra)^^j  - (ra)^lj>  4 O(C^)  (15c) 

The  first  term  on  the  right  side  of  Eqs.  (15)  is  the  contribution  of  the  pumping 

reaction,  and  the  second  term  is  the  modification  because  of  radiative  and 

collisional  deactivation.  The  pumping  reaction  increases  and  by  an 

amount  proportional  to  Y,  whereas  deactivation  increases  (or  decreases) 

and  by  an  amount  proportional  to  C Y.  Equations  (15)  are  similar  in  form 

4 5 

to  the  results  of  Holland  and  Mirels,  ’ where  a laminar  flame-sheet  solution 
was  obtained  by  expansion  about  x = o. 

Zero-power  number  density  and  gain  distributions  for  a cold-reaction 
HF  laser,  in  the  limit  Kj  — ► »,  Y = (C/Cq^*^*  given  in  Fig.  2.  For 


i 
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small  C,  and  vary  as  Y.  In  the  limit  C — ►«.  ^ Y,  and  for  v / 0, 

~ C The  population  inversion  is  total  for  small  C and  becomes  partial 

with  increases  in  J . The  gain  per  semichannel  is  a maximum  at  values  of 
C of  order  1.  In  interpreting  Figs.  2b  and  2d,  it  should  be  recalled  that,  in 
physical  variables,  the  gain  per  semichannel  is  proportional  to  a ~ (v  + 
for  a fixed  value  of  J.  Hence,  for  laminar  flow  and  fixed  J,  the  peak  gain  on 
the  2-1  transition  is  somewhat  larger  than  the  peak  gain  on  the  corresponding 
1-0  transition  (Fig.  2b),  although  this  is  not  true  for  turbulent  flow  (Fig.  2d). 

The  present  results  are  expected  to  be  modified  somewhat  by  vibration- 
vibration  collisions,  which  are  neglected  herein.  These  results  are  useful 
however,  for  providing  convenient  first  estimates  of  zero-power  number  density 
and  gain. 

C.  OSCILLATOR  CASES 

We  now  estimate  net  oscillator  output  power  following  a procedure 

7 8 

similar  to  that  used  by  Broadwell  and  Emanuel  and  Whittier.  A Fabry- 
Perot  optical  cavity  is  assumed.  Steady- state  lasing  is  assumed  to  be  initiated 

simultaneously  on  all  transitions  (v  = 0,  1, 2 Vj.-1)  at  some  station  denoted 

The  value  of  J is  assumed  to  be  the  same  for  each  lasing  transition,  i.e., 
A^  = A = const.  The  choice  for  J is  discussed  later.  All  photons  are  assumed 
to  have  the  s cune  ene  rgy,  i.e.,  6 ^ = f . 

8 

The  condition  for  steady-state  lasing  requires  that,  in  the  lasing  region. 


= [(-1)  In  (RjR2)]/[2n^^n^w]  v < v^ 


a a G 
o c 


(16) 
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is  the  number  of 


where  Rj  and  are  the  cavity  mirror  reflectivities,  and  n 


sc 


si'michanni'ls  bclwi'i'n  Ihc  mirrors.  'I’lio  ciiianlity  G o G /a  denotes  the 

• ' e V V o 


constant  value  of  G^  in  the  lasing  region.  Equation  (16)  does  not  apply  at 
V = Vj.  because  there  is  no  lasing  from  v^  + 1 — ► Vp  For  a given  optical  cavity, 
laser  geometry,  and  flow  conditions,  the  product  a G is  a constant.  However, 


o c 


the  individual  values  of  a and  G depend  on  the  choice  for  J.  Saturated  laser 

o c 


operation  corresponds  to  G — ► 0. 


llpstream  of  c;  .,  0 0.  The  variation  of  0 with  C.  downstream  of  f .,  is 

found  as  follows.  Equation  (5b)  can  be  expressed  as 


1 G = dN  /dC  - (dN  /dC)„  + rN  - r^.N^,+I  ,G  , 
vv  v"  v'^p  vv  v+1  v+1  v-1  v-1 


(17) 


Successive  solutions  of  Eq.  (17)  for  I^G^,  IjGj,  etc.,  yield 


[dNi/dC  - (dm/dOpl  - 


(18) 


The  rate  of  change  of  net  photon  output  is,  then,  after  algebraic  simplification 


Vf-1 


d0/dC=  E 1 G 

OV  V 


v = 0 


- clN  /clC  ) - r N 

V ^ V V ' 


(19) 
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Equation  (19)  can  be  deduced  more  directly  by  consideration  of  conservation  of 
vibrational  energy;  i.e.,  the  rate  of  change  of  output  power  equals  the  rate  of 
increase  of  vibrational  energy  created  by  the  pumping  process  less  the  rate  of 
increase  of  vibrational  energy  of  the  particles,  less  the  rate  of  energy  lost 
by  vibration-translation  collisions.  It  is  clear  that  vibration-vibration  collisions 
conserve  vibrational  energy,  and  if  these  had  been  included  in  the  present  study 
they  would  not  appear  in  Eq.  (19). 

We  now  express  N in  terms  of  G and  N_  which  are  known 
^ V o ToT 

quantities.  Successive  solutions  of  N = G , + AN  , for  N,,  N.,,  etc., 

vv-1  v-l  12 

indicate 


v-l 

N = ^ (G. /aM  + a'^N  (20) 

V . _ 1 o ' ' 

1-0 


where  the  summation  is  taken  equal  to  zero  for  v = 0.  An  expression  for  is 
found  by  summing  Eq.  (20)  for  v = 0,  1,  ...  v^.  The  result  is 


N 

o 


’f 

E a'' 


= N 


ToT 


E 

v = 0 


G 

V 


v^-  1 -v 

E A^ 


i=0 


(21) 


Substitution  of  Eq.  (21)  into  Eq.  (20)  yields 


N 

V 


= F N 

V 


ToT 


+ H G 

V c 


(22) 
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where 


= a''/  i a'' 

v = 0 


(23a) 


v^-1 


v^-l-v 


. v-1  a /a. 

= E (a^/a^)  z A‘ 

i=0  a'  ^ v=0  ° 1 = 0 


(23b) 


Equation  (23a)  can  be  expressed  as 


Vj.+  1 

- 1)/(A  - 1)  am 


= (v^  + 1) 


-1 


A = 1 


(24a) 


For  Eq.  (23b)  becomes 


v.+  l 


= [A^(l  + v^)/(A  ^ - 1)]  . (A  - 1) 


-1 


= V - (v^/2) 


A ^ 1 
A = 1 


(24b) 


f 

The  quantities  F and  H are  independent  of  C • Also,  E F = 1 E H = 0 

V V \r  * ^ \r  * 


v=0 


Thus,  is  known.  Substitution  of  Eq.  (22)  into  Eq.  (19)  yields 


v=0 


d0/dc  = C,dN^_,^/dC  - - 0,0^ 


(25a| 


where  it  is  assumed  that  dG^/dC  = 0. 
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r 


where 


C,  = i:  v(a  - F ) 

1 . V V 

v=  1 


f 


C,  = X r F 
2 V V 

v = l 


£ 


C-  = X r H 

3 V V 

v = 1 


If 


terms  of  order  A compared  with  one  are  neglected,  C,  = X va  , C,  = A, 

Vr  1 j V ^ 


Vf 


C-  = X r (o  /o  .).  The  latter  equals  v,  when  r = cj  , /a  = v.  From  the 
3 j V o v-1  ^ f V v-1  o 


definition  of  Cj,  it  is  clear  that  equals  the  number  of  photons  liberated  per 
lasing  species  molecule  (formed  by  the  chemical  reaction)  for  the  case  of  a 
saturated  laser  with  no  collisional  deactivation.  Similarly,  is  the  ratio  of 
the  net  collisional  deactivation  rate  in  a saturated  laser  to  the  collisional 
deactivation  rate  that  would  exist  if  all  the  lasing  species  molecules  were  in 


the  first  vibrational  level.  The  quantity  C_G  /N.^,  „ is  the  correction  to  C^, 

3c  ToT  Z 


i.e.,  increased  collisionad  deactivation,  due  to  a lack  of  optical  saturation. 

Equation  (25a)  applies  downstr  \m  of  Q The  net  photon  output,  up  to 
station  Q , is 


0 = C 


- C, 


C.  "2 


/ 


c 


^ToT^^ 


C3Gc(C 


-^i) 


(25b) 


■21 


1 


w 


The  optical  cavity  is  assumed  to  end  at  the  streamwise  station,  where  dp/ d^ 
becomes  zero  or  negative.  The  latter  station  is  denoted  C and  the  corre- 
sponding net  photon  output  and  power  output  are  denoted  and  P^,  respectively. 
If  AH  represents  the  chemical  energy  released,  per  mole  of  n^  particles,  the 
chemical  efficiency  of  the  laser  is 

n = P /(n  uwAH)  = (e/AH)p  (26) 

c e r e 


where  £ /AH  =1/3  for  a cold-reaction  HF  laser  operating  at  2.706  pm.  The 
product  2.327  p^  yields  the  output  at  this  wavelength  in  terms  of  kilojoules  per 
grzim  of  fluorine. 

Oscillator  performance  can  now  be  deduced.  Consider  the  limit 
Kj  — *■  with  a flame  sheet  given  by  Y ={C  /C  C < C Y = 1 for 

C > C Q.  Equations  (25)  yield 

C^(dp/da  = (Ci<C<Cj5)  (27a) 


^2^0  I C 
(m+l)C^  I C 


c =c 


(27b) 
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where,  for  (C^C  /mC^)  > [l  + (C  G /C 


c,C 


vm 


2^  e' 


mC , 


1 + 


^Z^dI^3^c 


2. 


2C, 


1/2 


'C.C^\ 


1/2 


2VD 


C_G 
3 c 


2C, 


(m  = 1 /2) 


(28a) 


= 1 


^2^0  ^3*^0 


Cl  C, 


(m  = 1)  (28b) 


.-1 


and  for  (C^C  ^^2^^ ' 


(28c) 


Equations  (28a)  and  (28b)  are  obtained  by  setting  Eq.  (27a)  equal  to  zero,  and 
correspond  to  cases  where  lasing  terminates  before  the  flame  sheet  reaches 
the  cliannel  centerline.  Equation  (28c)  corresponds  to  cases  where  lasing  is 
terminated  when  the  flame  sheet  reaches  the  channel  centerline.  The  quantity 


^1  can  be  found  from  Eq.  (16)  and  the  equations  given  in  Appendix  B.  In  most 


cases,  <<  1.  and  Eq.  (15c)  can  be  expressed 


/Cl' 

r,  r- 

a G 
o c 

[(«)v+l  - <”'v+2. 

- A 

(ra)  - (ra) 

L V v+lj 

Vdi 

"a  (a  , . - Aa  ) 

V V+1  V 

m + 1 

^v+1  ■ 

■ Aa 

V 

+ o(cf) 


(29) 
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1 


F 


whi'r*’  >•'*  nivt'ii  by  I'iq . (lb).  I0c|uali(>n  (Z'))  i.s  wrilt*'n  in  a form  (hat  can 

bo  solved  for  Q.  by  iforalion  with  C,.  0 used  as  the  inilial  esliinali'.  For  a 

Hiven  value  of  .1,  the  valui'  of  v that  yields  llu’  sniallesl  value  of  Qr.  i.e., 
corresponds  fc>  flu*  transition  with  tlu'  largest  /.I'ro  powt-r  (jain,  is  us«'d.  The 
choice  of  v = 1 (Z)  appears  appropriate  fi)i'  cold  - react  ion  lll''(ni'')  lasers.  If 
no  solution  for  f.  is  obtained,  tbreshohl  is  not  re;ubiMl. 

In  most  pr.ictical  di'vices,  the  optic. i|  medium  is  s.»t  i.  r.i  I I'd,  .ind  lasing 
i>rininate.s  near  (!  0-  lo  thes('  cases,  Feps.  (27)  an<i  (28)  can  lu'  considerably 

simplified.  Assume  that  (C./C, ..)*'=  0(G  )<<  1.  (C  Cl  /C' ,)''  <<  1 , and 
= d(l).  Eeptations  (27)  and  (28)  become,  for  (l,(;p/mClj  ^ 

It  1 (c,ci  /c:,)!"' 

IK.  4. 


(m  t 1) 


IntcJ  c:, 


t - m(m  It)- 

^mC.) 


n r'c 


TT 


1 C1((I  ) 

O 


I c)l(c:  ,ci^7c^)‘ 


(U)a) 


mC, 


1 - 


c 


c 


c 


and  for  1 - !•  < 


,-l 


( U)b) 


^'1 


, . /_L,  , , o(c-,  I 


m 1 I > I f- 


(Ulc) 
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(30d) 


Equations  (30)  are  plotted  in  Fig,  3 for  m = 1/2  and  define  saturated  laser 
performance. 

Values  of  Cj,  C^,  and  are  given  in  Table  1 as  a function  of  A for  values 
of  a , r^,  a^/a^,  and  v^.  appropriate  for  a cold- reaction  HF  laser.  For  this 
case. 


Cj  = 2.  11  - 


(31a) 


C2  = A(1  - A)(l  + 2A  + 3A^)/(1  - a"^) 


(31b) 


and,  to  within  1%,  = 3.0  + 1.7  A(1  - A).  It  is  seen  that  — ► 0 as  A — ► 0. 

This  result  is  explained  as  follows.  Recall  that  is  the  ratio  of  the  collisional 

deactivation  rate  in  a saturated  laser  to  the  collisional  deactivation  rate  that 

would  exist  if  all  the  lasing  species  molecules,  i.e.,  HF,  were  in  the  first 

vibrational  level.  When  A = 0 in  a saturated  laser,  stimulated  emission 

reduces  all  lasing  molecules  to  the  ground  state.  Hence,  there  is  no  collisional 

deactivation  in  this  case,  and  = 0.  This  limit  can  not  be  achieved  in  practice 

because  a — ► 0 as  A — ► 0 such  that  G = (a  G )/a  — ► « for  nonzero  values 
o coco 

of  a G . 
o c 

We  have  assumed  that  the  value  of  J is  the  same  for  all  lasing  transitions, 

but  the  choice  for  J has  not  yet  been  established.  In  Fig.  4 is  indicated  the 

effect  of  J on  HF  laser  output  for  the  realistic  range  of  conditions  0 < < 

5 2 

10  cm  /mole,  0.0  S < 20,  and  T = 300  K and  600  K.  For  these  conditions 
and  small  values  of  J,  the  optical  medium  is  saturated,  and  laser  output 


-25- 


Table  1.  Parameters  for  Defining  Oscillator  Performance 
for  Cold-Reaction  HF  Lasers 


A 

T 

^1 

^3 

0.00 

GO 

2.  1 10 

0.000 

3.  000 

0.  10 

1.151 

1.999 

1.  107  (-1) 

3.  159 

0.20 

8.  047  (-  1) 

1 . 866 

2.436  (-1) 

3.  291 

0.  30 

6.020  (-1) 

1.714 

3.959  (-1) 

3.  380 

0.40 

4.  581  (-  1) 

1 . 548 

5.616  (-1) 

3.424 

0. 50 

3.466  (-1) 

1 . 377 

7.  333  (-  1) 

3.422 

0.60 

2.  554  (-1) 

1 . 206 

9. 044  (-  1) 

3.  382 

0.70 

1. 783  (-1) 

1.041 

1.069 

3.  312 

0.80 

1.  116  (-  1) 

8.851  (-1) 

1. 225 

3.  220 

0.90 

5,  268  (-2) 

7.413  (-1) 

1.  369 

3.  1 14 

1. 00 

0.000 

6.  100  (-1) 

1 . 500 

3.  000 

%-  v;  Vj  - 3;  o Jo ^ = 1 + v;  a^,  a^,  a^,  a^  = 0.00,  0.17,  0.55,  0.28. 
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(bI 

T“ 

I — ^ — 1 — 1 — 1 — ^ — 
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/ 
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J 

^b 

^C„=10 
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Figure  4.  Effect  of  J Variation  on  Output  Power  from  Cold- 
Reaction  HF  Chemical  Laser.  Kj  m = 1/2, 

rv  = V,  = V + 1,  a^  = 0. 0,  0.  17,  0.  55,  0.28, 

OoGc  = cm'^/mole.  (a)  T = 300  K,  5 < < 20; 

(b)  T = 300  K,  0 < < 1;  (c)  T = 600  K,  5 < ; n 

<20;  (d)  T = 600  K,  ^ > Cr^Ki. 


increases  with  increase  in  J.  With  further  increase  in  J,  the  medium  becomes 
unsaturated  because  of  the  decrease  in  o^  with  increase  in  J,  i.e.,  Eqs.  (2c)  and 
(2f),  and  laser  output  reaches  a maximum.^  Thereafter,  an  increase  in  J results 
in  reduced  laser  output.  Further  increases  in  J result  in  threshold  not  being 

g 

reached.  Emanuel  and  Whittier  recommended  that  the  value  of  J that  yields  the 
maximum  output  power  be  used  to  estimate  laser  performance.  This  approach 
may  tend  to  overestimate  J and  laser  output  power  because  lasing  at  each 
streamwise  station  actually  occurs  on  those  transitions  with  the  highest  gain  — 
and  these  transitions  appear  to  be  at  somewhat  lower  values  of  J than  those  that 

g 

yield  the  highest  output  power.  The  method  of  Emanuel  and  Whittier  thus 

provides  an  upper  bound  on  laser  performance.  An  alternative  approach  is  to 

use  values  of  J that  are  observed  experimentally.  In  typical  cold -reaction  HF 

and  DF  lasers,  the  output  power  appears  to  be  centered  about  J = 5 to  8 and 

9 10 

J = 8 to  10  transitions,  respectively,  ' although  these  results  depend 
somewhat  on  initial  gas  temperatures  and  the  degree  of  diUition.  Hence,  the 
values  J = 7,  9 appear  to  be  representative  experimental  values  for  HF  and 
DF  lasers,  respectively.  Further  study  concerning  the  choice  for  J is 
desirable. 

J 


t 


The  maxima  in  Fig.  4 are  predicted  to  within  a few  percent  by  using  Cj  ” ® 
in  Eq.  (27  b). 
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m.  DISCUSSION  OF  RESULTS 


The  dependent  variables  are  functions  of  the  independent  variable  q; 
the  parameters  K^,  G^,  A,  J (or  T),  and  m;  and  lasing  molecule  distri- 
butive properties  (r^,  a^).  The  quantities  and  are  the 

ratios  of  streamwise  distance,  characteristic  diffusion  distance,  and  char- 
acteristic pumping  distance  to  characteristic  collisional  deactivation  distance, 
whereas  the  quantity  is  a measure  of  optical  saturation.  For  zero -power 
cases  with  j the  dependent  variables  N^,  G^  are  independent  of 

Cq.  In  oscillator  cases,  the  quantities  0g/Cj  and 

^1'  ^3  ^c^ ^2’  ^1^1^ ^2’  ^2  ^i^ cf  « > 

[l  + (C^  G^/C2)]  t the  quantity  (C^  ^ function  only  of 

In  practical  lasers,  the  assumptions  incorporated  in  Eq.  (30)  are 

approximately  valid.  The  resulting  variation  of  efficiency  and  output  power 

1 /2 

with  G2  Cq/iti  Gj  and  (G2  (G^  G^/C^)  is  shown  in  Fig.  3 for 

1 /2 

m = 1/2.  Note  that  (G2  (G^  G^/G2)  is  independent  of  pw.  For 

small  values  of  the  flow  corresponds  to  that  for  a premixed  laser.  In 
Fig.  3a,  it  is  indicated  that  laser  efficiency  decreases  as  increases.  The 
ordinate  in  Fig.  3b  is  proportional  to  output  power  per  semichannel  when  plenum 
temperatures  and  a stoichiometry  (p„/p)  ^.re  kept  constant.  It  is  shown  that 

r 

output  power  increases  with  until  G2  = [1  + (C^  G^/G2)]“*  and 

remains  constant  with  further  increases  in  Cj^.  In  the  latter  region,  lasing 
is  terminated  before  the  fuel  (H2)  reaches  the  centerline  of  the  oxidizer  (F) 
channel. 


il 

The  present  results  can  readily  be  expressed  in  terms  of  physical 
variables.  We  first  consider  a cold-reaction  HF  laser  with  combustor- 
generated DF  present  in  the  oxidizer  flow.  Let  p„  and  Pt^tt  denote  the 
initial  partial  pressures  of  atomic  fluorine  and  DF,  respectively,  and 
let  pp  characterize  the  mean  value  of  p^^p  in  the  reaction  zone.  We 
assume  that  HF  and  DF  are  the  major  deactivators  of  HF(v)  and  use  only 

the  leading  terms  in  Eqs.  (A- 18)  and  {A-19).  Recall  that  C = k^Vx/u. 

cd 

It  follows  that 


I' 


f 


0.0752  p(torr)^(cm>  l±J<^\  [i  , 0.6  {Pp^/p„)l 
(p/Pp)(T/400)^  \ / DF  F 


(32) 


L 1/2 

Characterize  laminar  and  turbulent  flame  shapes  by  y = A^  (Sl.c/u)  ' and 

D D 

T 3/2 

y^  = A.pX,  respectively,  where  Aj^  = 0(1),  A.^,  = 0(0.  1),  and  2)  = CT  /p  is 

L 2 T 

the  laminar  diffusion  coefficient.  Then,  Xj^  = (w/Aj^)  n/U>,  Xj^  = w/A.^,,  and 
the  corresponding  laminar  and  turbulent  values  of  are 


. L.1/2  _ 2.072  p(torr)w(cm)  (z.  88  X lO"^  \ / 2 \ r,  . . , , 1 


/2 


(33a) 


-4 


where  C = 2.88  v 10  corresponds  to  the  diffusion  of  into  He  and 


jT  ^ 0._752,ELt.orrMaT,>  ^[,+0.6  (p„,,/p,,))  (33b| 

° (p/Pp)(T/400)^  u A,^  DF  F 
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Similarly,  if  we  let  c = a . and  n = p„/RT,  Eq. 

O O,  J T I? 

lasing  molecule. 


0.234  (p/pj,)(T/400)^^^  In  (RiR2)'^ 
” n^^  p(torr)w(cm) 


(16)  becomes,  for  a HF 


exp[J(J-l)Tj^/T] 
J(1  + 0.063  J) 


(33c) 


where  a - has  been  evaluated  from  Eqs.  (2).  Note  that  G is  independent 
of  pw.  Recall  that  -2.11,  = 3,  and  = exp  (-2  JT^^/T)  for  a HF  laser 

with  Hence,  the  dependence  of  laser  performance  on  physical  vari- 

ables can  be  readily  deduced,  c.  g.  , Fig.  3 and  Eqs.  (33). 

We  now  note  explicit  expressions  for  <p^  and  for  a saturated  laser 
(G^  = 0)  in  the  limit  Kj — ► C2Cj^/mCj>l.  From  Eqs.  (30)  and  (33) 

(p/p_)^/^(T/400)’^/^(C/2.88  x 10’^)^^^  (A.  /2) 

0 = 0.  228  -fry £ . (34a) 

p(torr)w(cm)[l  + (0. 6 Pj^p/Pp)]  ' 


when  m = 1 /2  and 


0 


e 


1 

0.665 

^2 


(p/Pp)(T/400)^(u/4  X 10^)(A.j,/0.  1) 
p(torr)w(cm)ll  + (0.6  Pj^p/pp)J 


(34b) 


when  m = 1.  The  dependence  of  on  physical  variables  is  considerably 
different  for  the  laminar  and  turbulent  mixing  cases  except  for  the  dependence 
on  pw.  The  corresponding  length  of  the  lasing  zone  is,  from  Eqs.  (30b)  and 
(32) 
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, , me,  (p/pp)(T/400)2  (u/4xl0^ 

- 0-2  0.0752  p(torr)  [1  + (0. 6 

The  lasing  length  is  proportional  to  u/p,  which  is  expected  from  binary  scaling 

considerations,  and  is  independent  of  w. 

Similar  expressions  can  be  deduced  for  a cold-reaction  DF  laser  with 

combustor-generated  HF  present  in  the  oxidizer  flow.  For  a DF  lasing 

molecule,  a = 1,  T„  = 15.6  K,  M = 21  g/mole,  and 
K. 

, X 10^^  = 0.662  [l  + 0.045  j]  erg-cm^  (36a) 

O y J 

j/Mq  j = 1 + V - 0.005  (36b) 

a . 2.03  X 10^  (400/T)^''^/mole  (36c) 

Equations  (36)  are  correlations  of  estimates  by  Herbelin  that  are  believed  to 

be  accurate  to  within  about  10%  for  1 ^ J S 13,  v < 6.  Consider  HF  and  DF 

to  be  the  major  deactivators  of  DF(v),  approximate  Pj^p  by  the  initial  value 

of  Pp,  and  include  only  the  leading  terms  in  Eqs.  (A-19)  and  (A-20).  Equations 

(32)  through  (35)  then  apply  for  the  cold  - reaction  DF  laser  if  [ 1 +0.6  (p  /p  )] 

is  replaced  by  0. 4 [ 1 + 1.6  (pj^p/Pp)],  if  [1  + 0.063J]  is  replaced  by  0.365 

-4 

[1  + 0.045  j],  and  if  C = 2.  35  X 10  (corresponding  to  diffusion  of  D2  into  He) 
is  used  in  Eq.  (34a).  For  a cold-reaction  DF  laser,  a^  = 0,  0.  10,  0.24, 

I 0. 38,  and  0.  28;  hence,  = 2.  84  and  = 4 for  C2  = exp  (-2  JTj^/T)  1. 

The  product  1.  641  0^  yields  the  laser  output  in  terms  of  Icilojoules  per  gram 
of  fluorine  for  a mean  output  wave  length  of  3. 837  pm. 

i 

I 


(35) 
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The  two-level  model  of  Mirels,  Hofland,  and  King  corresponds  to  the 

choice  V.  = a =r.=A=a  /a.  = l.  For  this  case,  C,  = C_  = C_  = 1/2,  and 
X Z X X ^ X ^ 

the  present  results  agree  with  those  of  Ref.  6 For  a given  lasing  species, 
scaling  laws  deduced  from  the  multiple-level  model  differ  from  those  deduced 
from  the  two-level  model  only  in  the  dependence  of  C^,  and  on  T and 

J.  The  effect  of  the  latter  on  laser  performance  is  approximated,  in  the 
two-level  model,  by  appropriate  normalizations;  e.  g,  , chemical  efficiency 
for  Cj-j  / 0 is  normalized  to  the  value  for  = 0,  and  by  the  use  of  suitable 
mean  rates;  e.g.  , was  used  to  characterize  HF  collisional  deactivation. 

7 

The  mixing  model  used  by  Broadwell  is  illustrated  in  Fig.  5.  It  is 
shown  in  Appendix  C that  the  latter  leads  to  results  that  are  identical  with 
the  flame -sheet-mixing  model  used  herein  if  mean  chemical  rate  are  used 
in  both  models. 


^The  present  oscillator  solution  becomes  identical  with  the  solution  of  Ref.  6 
if  the  present  variables  C/Cj,  ^3  ^c^^2’ 

identified  with  the  variables  Q,  T],  Gg/<7w[F]^,  and  Kj  in  Ref.  6. 


k. 


-35- 


IV.  CONCLUDING  REMARKS 


I 

i 

i 


It  is  clear  from  Eqs.  (34)  that  the  dependence  of  0^  on  physical 

variables  is  different  for  laminar  and  turbulent  flows.  Hence,  it  is  important 

to  establish  whether  the  flow  is  laminar  or  turbulent  before  applying  the 

prjpsent  results.  The  Reynolds  number  per  unit  length  at  the  oxidizer  nozzle 

11  3 

exit,  for  a typical  current  laser  design  is  of  order  R^/cm  = 0(10  ),  e.g.  , 

5 

helium  at  p = 10  torr,  T = 300  K,  u = 4 x 10  . Most  of  the  energy  is  ex- 
tracted within  a few  centimeters  of  the  nozzle  exit.  This  suggests  that 

the  mixing  is  laminar  in  the  region  of  interest,  as  was,  in  fact,  observed 
12  13 

by  Varwig  and  by  Shackleford  et  al.  Shock-induced  separation  at  the 

nozzle  exit  and  blunt  nozzle  trailing  edges  may  hasten  transition  to 
12 

turbulence.  In  the  latter  case,  the  mixing  in  the  lasing  region  may  be 
classified  as  "transitional.  " In  these  cases,  the  choice  m = 1/2  or  an 
alternative  choice  for  y^  = yj(x)  may  be  more  appropriate  than  the 
choice  m = 1,  which  implies  fully  developed  turbulent  flow.^  Further  study 
is  needed. 

It  should  also  be  noted  that  the  simplified  models  of  Mirels,  Hofland, 

6 7 

and  King,  Broadwell,  and  the  present  study  neglect  features  of  the  flow 


Note, from  Eq.  (35), that  Reynolds  number  based  on  lasing  length  is  independent 
of  pressure  level.  Hence,  an  increase  in  pressure  level  does  not  necessarily 
result  in  turbulent  mixing  in  the  lasing  region. 
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that  may  be  important  in  some  cases.  In  particular,  initial  flow  nonuniformity 
(because  of  nozzle  wall  boundary  layer);  pressure  nonuniformity,  i.  e.  , 
shocks,  edge  effects,  and  lateral  gas  expansion;  and  heat  addition  effects, 
i.e.  , streamwise  temperature  variation  effect  on  rates,  are  neglected. 
Appropriate  mean  values  are  needed. 
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APPENDIX  A.  REACTION  MODEL 


1 

The  simplified  reaction  model  used  in  this  report  is  illustrated  by 
considering  pumping  and  colHsional  deactivation  reactions  for  a HF  laser. 

Backward  reactions  are  neglected. 

A.  PUMPING  REACTION 
1.  COLD -REACTION  LASER 

The  cold-reaction  HF  laser  uses  molecular  hydrogen  and  atomic  fluo- 
rine. Vibrationally  excited  HF  is  generated  by  the  pumping  reaction 


F t H^  -14.  HF(v)  + H 


(Al) 


(v)  . — f "771 

where  k^  is  the  rate  coefficient  for  production  of  HF(v),  and  ~ L k'  ^ 

is  the  overall  rate  coefficient.  The  rate  of  change  of  n^,  along  a single 

stream  tube  with  velocity  u,  assuming  a hydrogen  rich  mixture,  is 


u d n„/dx  = - k n„  n_ 
F p H , F 


3 - k n_ 
P F 


(A2a) 

(A2b) 


-1 


The  quantity  k s k n„  has  the  units  sec  , and  a mean  value  is  used.  The 
P P “2 

reciprocal  of  kp  is  the  characteristic  pumping  time.  If  the  reaction  starts  at 
X with  an  initial  concentration  (n_)  (Fig.  1),  the  downstream  concentration 

O f OD 

of  atomic  fluorine  along  the  stream  tube  is 


np/(nF)^  = e 


-kp(x-Xo)/u 


(A3) 
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Let  Vk  and  note  X]  ^ From  Eqs,  (Al)  and  (A2),  the 

vibrationally  excited  HF  created  by  the  pumping  reaction  is,  for  each 
stream  tube, 

u(d  n /dx)  = a k n^-,  (j 

V p V p r ' 

The  net  rate  of  creation  of  HF(v)  at  a given  streamwise  station,  due  to  all 
the  stream  tubes,  is  found  by  integrating  Eq,  (A4)  between  y = 0 and  y£(x). 
The  result,  after  introduction  of  nondimensional  variables,  is 

(d  N^/dt)p  = % Kl  Nf 

where 

K,  Nf.  = Kj  e-KlC  (dY^/dc^)  dj^  U 

= dY/dc  (K^--) 

and  n^  = (n^)^,  = y£(x^)/w,  Y = y^Cxl/w. 

2.  CHAIN-REACTION  LASER 

The  chain- reaction  HF  laser  uses  H2  and  partially  dissociated  F^. 
Vibrationally  excited  HF  is  created  by  the  chain  reaction 

k '"f 

H + F -Pi^  E [HF(v)l  + H (A 

v=0 

k . '"f 

H + F2  E [HF(v)]£^  + F (A 


where  subscripts  c and  h denote  cold  and  hot  reactions,  respectively,  and 

V£ 

k = Y , and  k . = [Subscript  c was  omitted  in  Eq.  (Al).] 

P'^  v=0  v=0  P»" 

The  rate  of  change  of  species  along  a stream  tube  is  found  from 


^2  ^Tot\  _ *^h\ 

dx  \ dx  I ‘V  dx  ) -\  dx  I 


n,-  •^TT 

p.c  F H; 


d n. 


dx 


"Tot\  ''f\ 

\ k~'  \ h 


k , n„  n„ 
p,h  H F2 


(A7a) 


(A7b) 


where  n_  ™ n . It  follows  from  Eq.  (A7) 

ToT  _ V T ' ' 

v=0 


“f  ^ "h  = <"f 


(A8) 


The  solution  of  Eqs.  (A7)  is  found  by  using  an  approach  similar  to  that  used 
14 

by  Skifstad  for  the  premixed  laser  case.  In  order  to  obtain  a closed-form 
solution,  we  assume  that 


k 

P.  c 


”f 


) 


It  follows  that,  along  each  streamtube 


(A9a) 

(A9b) 


(AlOa) 

(AlOb) 
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Define  ^ ^ The  production  of  n^  by  the 

pumping  reaction  is,  then. 


u(d  n /dx)  = a k n„ 

V p V p 


(All) 


where  a - a + a , . The  net  production  of  n at  each  streamwise  station 

is  found  by  integrating  Eq.  (All)  between  y = 0 and  y^(x).  Introduction  of  non- 

dimensional  variables,  with  n = (n_  ) , yields 

r ^2’^ 


(d  N^/dC)p  =a^K^ 


(A12a) 


where 


•'l  V = *^1  /'  dc^ 


(A12b) 


Equations  (A12)  have  the  same  form  as  Eqs.  (AS).  However,  a , n . and  k 

V r p 

are  defined  differently. 

The  present  solution  is  self  consistent  provided  (n_)  /(n,,)  = k , / 

r 00  H oo  p»  h 

kp,c'  most  HF  lasers,  (n^)^  ^ 0,  (nj^)^  = 0.  In  these  cases,  the  present 
solution  becomes  applicable,  after  a short  transient,  for  flows  with  small 
initial  F^  dissociation.  Values  of  n„  and  n„,  after  the  transient,  are 


and  the  pumping  rate  for  the  chain  reaction  is,  then. 


h (n„) 
k = p,  c ' F <» 

P [1  + (k  /k  . )] 

^ p,c  p,h'' 


(A14) 


Denote  the  initial  dissociation  level  by  a = [n_/(n_  + 2n„  )]  . The  ratio  of 

^ jc  i?  * 2 * 

the  pumping  rate  in  a chain-reaction  laser  to  the  rate  in  an  equivalent  (same 
net  fluorine  flow)  cold-reaction  laser  is  of  order  a-  For  low  values  of  the 
parameter  may  be  reduced  to  values  of  order  one,  and  the  chain-reaction 
laser  will  have  considerably  lower  chemical  efficiency  than  an  equivalent 
cold-reaction  laser. 

B.  COLLISIONAL  DEACTIVATION 

We  now  consider  vibration-translation  (v-t)  deactivation  of  excited 
species  and  neglect  vibration -vibration  transfer.  The  v-t  deactivation  of 


n^  by  species  n^,  namely. 


is  found  from 


n + n.  ► n , + n. 

V 1 v-1  1 


u d n^/dx  = - n.  n^ 


(A15) 


The  deactivation  of  n^  by  all  species  can  be  expressed 


u d n /dx  = - k' j'  n 
V cd  V 


(Alba) 
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where 


^cd  ^ 

and  = 0.  The  quantity  is  assumed  to  be  a constant,  and  characteris- 
tic values  of  the  major  deactivators  are  used  in  Eq.  (Al6b).  For  arc-driven 
cold-reaction  HF  lasers, the  major  deactivator  is  HF  (and  possibly  H,  F). 
The  number  density  of  each  of  the  latter  can  be  characterized  by  (n_)  . For 

P 00 

combustion -driven  HF  lasers  the  introduction  of  major  deactivation  species 
by  the  combustor,  e.g.,  DF,  needs  to  be  considered. 

C.  RATE  COEFFICIENTS 

Pumping  and  collisional  deactivation  rates  for  a HF  laser,  from 

15  —7, 

Cohen,  are  repeated  here  for  convenience.  Units  are  T =®K  and  k = cm  / 
mole-sec.  The  overall  pumping  reaction  rates  [Eq.  (A6)]  are 


ir  =1.6X10‘^e-8<>5/T 

p,  c 

k , = 1.  2 X '1208/T 

p,h 

and  the  distributional  coefficients  are 


(A17a) 

(A17b) 


av,c  = 0,0.  17,0.55,0.28  (v  = 0,1, 2, 3)  (A17c) 

®v,h  ^ 0,0,0,0.08,0.  13,0.35,0.44  (v  = 0,  1 , . . . , 6)  (A17d) 


The  major  collisional  deactivation  process  in  arc -driven  lasers  is 


HF(v)  + HF  - HF(v  - 1)  + HF 
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for  which  the  rate  coefficient  is 


14-1  42 

10  T + 3.  5 X 10  T 


. 26j 


(A18) 


The  latter  differs  somewhat  from  the  rate  used  by  Mirels,  Hofland,  and 

King.^  Collisional  deactivation  of  HF(v)  by  DF,  in  combustion -driven  lasers, 

9 

can  be  estimated  from 


= v(1.9  X 10^"^  + 1.35  X 10^  T^] 


(A19) 


Equation  (A17)  neglects  vibration-vibration  collisions,  which  may  be 
significant. 

The  rate  coefficient  for  the  collisional  deactivation  of  DF(v)  by  DF  is 


v|1.2X  10^^  T'^  + 1.4X  10^  T^ 


.96) 


(A20) 


The  rate  coefficient  for  the  deactivation  of  DF(v)  by  HF  is  the  same  as  that 
indicated  in  Eq.  (A  19). 


N.  Cohen,  The  Aerospace  Corporation,  private  communication. 


-45- 


.4 


1 


t; 

’•ti 

APPENDIX  B,  ZERO-POWER  NUMBER  DENSITY  DISTRIBUTION 

V- 

Number  density  distributions,  in  the  absence  of  radiation,  are  found  i 

herein  for  the  case  Y = (C/Cjj)”^,  - ®,  r^  = v,  and  <3.  Integration  of 


TURBULENT  FLOW  (m  = 1) 


^3  ~ ® (B2a) 

^2  ^ ^2  ■'■  ®3  ■ ^^2  ■'■  3a2)e"^^  + 2a^e~^^  (BZb) 

Cj-jN^  = a^  + a^  + a^  - (a^  + 2a2  + 3a^)  e"^ 

+ (^2  ■*■  ^^3^  (B2c) 

^0^0=^  ■ + ^2  + ^3)  (B2d) 
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APPENDIX  C.  MIXING  MODELS 


» 


It  is  now  shown  that,  when  mean  values  of  reaction  rates  are  used,  the 

7 

scheduled-mixing  model  of  Broadwell  is  equivalent  to  the  mixing  model  of 
Mirels,  Hofland,  and  King^  used  herein. 

The  mixing  model  of  Broadwell  is  illustrated  in  Fig,  5.  For  a HF 
laser,  the  quantities  w^  and  w^  represent  the  semichannel  widths  of  the  hydrogen 
and  the  fluorine -diluent  streams,  respectively.  Both  streams  have  the  same 
velocity  u.  The  reaction  zone  is  assumed  to  be  bounded  by  two  similarly 
shaped  curves,  denoted  by  y^  and  y^  in  Fig.  5.  Fluid  properties  in  the  re- 
action zone  depend  only  on  x.  Free- stream  particles  that  intercept  the 
bounding  curves  between  station  x and  station  x + Ax  are  assumed  to  be  uni- 
formly distributed  throughout  the  reaction  zone.  The  axial  distance  required 
for  all  the  reactants  to  enter  the  reaction  zone  is  denoted  by  x^^  and  is 
equivalent  to  the  diffusion  distance  in  the  flame  sheet  model.  Note  that  w^ 
and  are  the  values  of  y^  and  at  x^.  Although  Broadwell  considered 
only  turbulent  mixing  with  w^  = w^,  we  only  require  that  y^  and  be  similar 
in  shape,  i.e.,  y ^^2.  ~ constant.  Introduce  change 

of  any  species  n.  in  the  reaction  zone  can  be  expressed  as 


dx  ^ i 2 dy,  ^ 1 dy,  ^ y,  dx  ” y,  dx 


dy^  n.u  dy^ 


+ (U: 


(Cl) 


where  (n.)j  and  (n^)2  represent  the  free-stream  values  of  n.  intercepting  curves 
y^  and  y^,  respectively.  The  first  term  on  the  right  side  of  Eq.  (Cl)  represents 


!l 
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the  rate  of  increase  of  n.  because  of  particles  entering  the  reaction  zone; 
the  second  term  represents  the  decrease  in  n^  because  of  the  lateral  expan- 
sion of  the  reaction  zone;  the  third  term  represents  the  rate  of  change  of 
n.  because  of  chemical  and  radiative  reactions.  Since  and  are  assumed 
to  be  similar  in  shape,  a mean  free- stream  value  of  m,  denoted  (n. ) can 
be  introduced  such  that 

(n.)^  = [(n.)2y2  - (ni)jyj]/yj  (C2) 

Equation  (Cl)  can  then  be  written  as 


u 


d([  n.  - (n.)^  ] y^} 
dx 


= 0). 
1 


(C3) 


For  lasing  molecules,  the  right  side  of  Eq.  (C3)  has  the  same  terms  as  the 

right  side  of  Eq.  (3).  Introduce  w = w-  - w,,  N = n y-/w  n , G = g y,/ 

Z 1 V v-'f  r V ‘’v^f 

awn,  and  note  (n  ) =0.  Equation  (C3)  then  reduces  to  Eqs.  (5).  The 

Vr  yOO  -I'/ 

term  (dN  /d  C)  in  Eqs.  (5)  is  found  by  integrating  Eq.  (C3)  for  n.  = n_.  The 

V p 1 F 

results  are  the  same  as  those  deduced  in  Appendix  A. 

, 7 

i is,  the  mixing  model  of  Broadwell  corresponds  to  a flame  sheet 
model  v-it’h  premixed  reactants  [ the  number  density  being  given  by  Eq.  (C2)]  , 
a flame  sheet  location  y£  “ y2  “ ^ diffusion  distance  Xj^,  which  is 

evaluated  by  consideration  of  the  actual  diffusion  process.  At  each  axial 
station,  number  densities  in  the  Broadwell  model  correspond  to  average  values 
in  the  flame- sheet  model.  The  equivalence  requires  that  the  coefficients  of 
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SYMBOLS 


A,  A , A . 

\'  V , vT 


Cj,  C2,  Cj 


Ky*  Kv,  J 


I . I T 

V V,  J 

I , I T 

V V,  J 


-ZJ  T../T 

= o K 

= friiction  of  txirticlos  pumpod  into  lovol  v,  k^'^Vk 

P P 

= functions  of  J Tj^/T,  Eq.  25 

= defines  in  satvirated  oscillator,  Eq.  (22) 

= normalized  gain  per  seniichannel  in  optical  cavity, 
o G /a 

V V o 

- normalized  gain  per  semichannel,  Eq.  (4) 

= gain  per  unit  length  on  transitition  v + 1,  J - 1 -•  v,  J,  cm' 
= defines  effect  on  N of  G / o,  Eq.  (22) 

= net  energy  (joules)  released  by  pumping  reaction  per 
mole  of  n 


= Local  lasing  intensity,  watt/cm 
= normalized  value  of  lasing  intensity,  Eq.  (4) 

= rotational  quantum  nuniber 

= ratio  of  net  pumping  rate  to  v = 1 -•  v = 0 deactivation  rate, 

k /k 
p cd 

- 1 

= mean  rate,  sec 

3 

= rate  coefficient,  cni  /(mole-sec) 

= molecular  weight  of  excited  species 

= 1/2,  1 for  laminar  and  turbulent  flames,  respectively 


ToT 


= measure  of  net  amount  of  particles  in  levels  v = 0 to  v = 
at  station  x 

= measure  of  net  amount  of  species  i at  station  x,  Eq.  (4) 

3 

n^  = local  number  density  of  species  i,  mole/cm 

n^  = reference  number  density  (equals  [np]^  and  [n  ] for  cold 

reaction  and  chain  reaction  HF  lasers,  respectively) 

n = number  of  semi  channels 
sc 

P = net  output  power,  per  semi  channel,  emitted  up  to  station  x, 
watt /cm 

= net  oscillator  output  power,  per  semichannel,  watt/cm 
R2  = reflectivity  of  mirrors  1 and  2 

r = relative  deactivation  rate, 

V ’ \ cd  sp/'  cd 

T = rotational  and  translational  temperature,  °K 
= characteristic  rotational  temperature,  °K 

u = flow  velocity  in  x direction,  cm/sec 

V = vibrational  level,  v = o,  1,  . . . , v^ 

v^  = highest  vibrational  level  considered 

w = oxidizer  channel  semiwidth,  cm  (Fig.  1) 

X = axial  distance,  cm 

Xj^  = characteristic  diffusion  distance,  y^(Xj^)  = w 

X -flame  sheet  abscis  sa,  cm 
o 


1 


Y = yj(x)/w 
Y,,  =yj(x„l/w 

yj(x)  = flame  sheet  ordinate,  cm 
e = mean  value  of  e 

V 

= energy  (joules)  per  mole  of  photons  for  transition 
v+1,  J - l-*v,  J 

0 = net  photon  output  up  to  station  Q,  per  initial  reference  (n^) 
Vf-1 

particle,  0 = ^ 0^ 

v=0 

0^  = value  of  0 for  single  transition  v+l,J-l-*v,  J 

0^  = net  photon  output  per  initial  n^  particle  for  an  oscillator 

C = normalized  axial  distance,  x/u 

= ratio  of  characteristic  diffusion  distance  to  characteristic 
deactivation  distance, 

= value  of  C at  end  of  lasing  zone  in  oscillator 

Q = value  of  Q at  which  lasing  is  initiated 
i 

2 

a , c , = cross  section  for  stimulated  emission,  Eq.  (2),  cm  /mole 

V J 

1)  = chemical  efficiency,  Eq.  (26) 

c 

Subscripts 

cd  = collisional  deactivation 
F,  F2  = atomic  or  molecular  fluorine,  respectively 
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r 


p 

sp 


V or  V,  J 


= species  1 
= pumping  reaction 
= spontaneous  emission 
= vibrational  energy  level 
= associated  with  transition  v + 
= free- stream  value 
= corresponding  to  v = 0 or  x 
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THEJVANA.  GETTING  LABORATORIES 

The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver* 
satility  and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory;  Launch  and  reentry  aerodynamics,  heat  trans- 
fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory;  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory;  Development  of  new  materials:  metal 
matrix  composites  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue-induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory;  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere:  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  systems. 
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THE  AEROSPACE  CORPORATION 
El  Segundo,  California 


